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We review the key features of POWHEG, a method for interfacing parton-shower Monte Carlo generators to fixed 
next-to-leading order QCD computations. We describe a recently introduced framework, the POWHEG BOX, that 
allows the automatic POWHEG implementation of any given NLO calculation. We present a few results for Higgs 
boson production via vector boson fusion and Z + 1 jet production, both processes available in the POWHEG BOX. 



1. The POWHEG method 

In the past two decades, next-to-leading or- 
der (NLO) QCD computations have become stan- 
dard tools for phenomenological studies at lep- 
ton and hadron colliders. QCD tests have been 
mainly performed by comparing NLO results with 
experimental measurements, with the latter cor- 
rected for detector effects. 

On the experimental side, leading order (LO) 
calculations, implemented in the context of gen- 
eral purpose Shower Monte Carlo (SMC) pro- 
grams, have been the main tools used in the anal- 
ysis. SMC programs include dominant QCD ef- 
fects at the leading logarithmic level, but do not 
enforce NLO accuracy. These programs were rou- 
tinely used to simulate background processes and 
signals in physics searches. When a precision 
measurement was needed, to be compared with a 
NLO calculation, one could not directly compare 
the experimental output with the SMC output, 
since the SMC does not have the required accu- 
racy. The SMC output was used, in this case, to 
correct the measurement for detector effects, and 
the corrected result was compared to the NLO 
calculation. 

In view of the positive experience with QCD 
tests at the NLO level, it has become clear that 
SMC programs should be improved, when possi- 
ble, with NLO results. In this way a large amount 
of the acquired knowledge on QCD corrections 
would be made directly available to the experi- 
mentalists in a flexible form that they could easily 
use for simulations. 



The problem of merging NLO calculations with 
parton shower simulations is basically that of 
avoiding overcounting, since the SMC programs 
do implement approximate NLO corrections al- 
ready. 

The features implemented in POWHEG can be 
summarized as follows: 

- infrared-safe observables have NLO accu- 
racy; 

- collinear emissions are summed at the 
leading-logarithmic level; 

- the double logarithmic region (i.e. soft and 
collinear gluon emission) is treated correctly 
if the SMC code used for showering has this 
capability. 

In the case of HERWIG J1I2J this last requirement 
is satisfied, owing to the fact that its shower is 
based upon an angular-ordered branching. 

In ref. j3] a method, called POWHEG (for Pos- 
itive Weight Hardest Emission Generator), was 
proposed that produces positive-weighted events, 
and that is not SMC specific. In the POWHEG 
method the hardest radiation is generated first us- 
ing the exact NLO matrix elements. The POWHEG 
output can then be interfaced to any SMC pro- 
gram that is cither p T -ordcrcd, or allows the im- 
plementation of a p T vetoo However, when inter- 
facing POWHEG to angular-ordered SMC programs, 



X A11 SMC programs compatible with the Les Houches In- 
terface for User Processes [1] should comply with this re- 
quirement. 



the double-log accuracy of the SMC is not suffi- 
cient to guarantee the double-log accuracy of the 
whole result. Some extra soft radiation (techni- 
cally called vetoed-truncated shower in ref. [3]) 
must also be included in order to recover double- 
log accuracy. In fact, angular ordered SMC pro- 
grams may generate soft radiation before gen- 
erating the radiation with the largest p T , while 
PQWHEG generates it first. When POWHEG is in- 
terfaced to shower programs that use transverse- 
momentum ordering, the double-log accuracy 
should be correctly retained if the SMC is double- 
log accurate. The ARIADNE program [5] and 
PYTHIA 6.4 [6J (when used with the new showering 
formalism), both adopt transverse-momentum or- 
dering, in the framework of dipole-shower algo- 
rithm, and aim to have accurate soft resummation 
approaches, at least in the limit of large number 
of colours. 

In the POWHEG formalism [TJ, the generation of 
the hardest emission is performed first, according 
to the distribution 



da = B($ B ) d$ B 
B(<P b ) 



A (p- in ) 

A (fcr ($r)) d$ rad 



(1) 



where B ($b) is the leading order Born contribu- 
tion and 



B(^b) = B(^ b )- 



V($ B )+ d$ iad R($ R ) 



(2) 



is the NLO differential cross section used to gener- 
ate the Born variables (V ($b) and R ($_r) stand 
respectively for the virtual and the real correc- 
tions), and 
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is the POWHEG Sudakov form factor. The trans- 
verse momentum of the emitted particle is here 
indicated with fc T ($_r). The cancellation of soft 
and collinear singularities is understood in the 
expression within the square bracket in eq. @. 
Partonic events with hardest emission generated 



according to eq. (p} are then showered with a fc T - 
veto on following emissions. For all the techni- 
calities concerning this and other issues, we refer 
to [3T7] . 

Up to now, the POWHEG method, in the con- 
text of hadron colliders, has been applied to ZZ 
pair hadroproduction [8], heavy-flavour produc- 
tion [9], Drell-Yan vector boson production p~0| 
111] , Higgs boson production via gluon fusion |12[ 
113] . Higgs boson production associated with a 
vector boson (Higgs-strahlung) [13 , single-top 
production [14] , Higgs boson production in vector 
boson fusion [15] and Z + 1 jet production [16]. 

2. The POWHEG BOX 

The POWHEG BOX [17] is a computer framework 
that implements in practice the theoretical con- 
struction of POWHEG. The aim of the POWHEG BOX 
is to construct a POWHEG implementation of a NLO 
process, given the following ingredients: 

1. The list of all flavour structures of the Born 
processes. 

2. The list of all flavour structures of the real 
processes. 

3. The Born phase space. 

4. The Born squared amplitudes B, the color 
correlated ones Bij and spin correlated ones 
Bfj, u . These are common ingredients of NLO 
calculations performed with a subtraction 
method. 

5. The real matrix elements squared for all rel- 
evant partonic processes. 

6. The finite part of the virtual corrections 
computed in dimensional regularization or 
in dimensional reduction. 

7. The Born colour structures in the limit of a 
large number of colours. 

With the exception of the virtual corrections, 
all these ingredients are nowadays easily ob- 
tained. A matrix element program can be used 
to obtain Qj and ([5]). The colour-correlated and 



spin-correlated Born amplitudes can also be gen- 
erated automatically. Recent progress in the au- 
tomatization of the virtual cross section calcula- 
tion may lead to developments where even the vir- 
tual contribution (j6j) may be obtained in a pain- 
less way. Given the ingredients listed above, the 
PQWHEG BOX does all the rest. It automatically 
finds all the singular regions, builds the soft and 
collinear counterterms and the soft and collinear 
remnants, and then generates the radiation using 
the POWHEG Sudakov form factor. 



3. Higgs boson production in vector-boson 
fusion 

Higgs boson production via vector-boson fu- 
sion (VBF) is expected to provide a copious 
source of Higgs bosons in pp-collisions at the 
Large Hadron Collider (LHC) at CERN. It can 
be visualized as the inelastic scattering of two 
quarks (antiquarks) , mediated by t-channel W or 
Z exchange, with the Higgs boson radiated off the 
weak bosons. It represents (after gluon fusion) 
the second most important production process for 
Higgs boson studies [18)19] . Once the Higgs bo- 
son has been found and its mass determined, the 
measurement of its couplings to gauge bosons and 
fermions will be of main interest [20 21 . Here 
VBF will play a central role since it will be ob- 
served in the H ^ tt [22125] . H ->• WW [23125] 
and H — > 77 [25] channels. This multitude of 
channels allows to probe the different Higgs boson 
couplings. In addition, in order to distinguish the 
VBF Higgs boson signal from backgrounds, strin- 
gent cuts are required on the Higgs boson decay 
products as well as on the two forward quark jets 
which are characteristic for VBF. The efficiency of 
these cuts has to be evaluated on the basis of the 
most updated simulation tools and experimental 
inputs. 

An additional feature characterizing VBF 
Higgs boson production is the fact that, at lead- 
ing order, no colored particle is exchanged in the 
t channel so that no i-channel gluon exchange is 
possible at NLO, once we neglect the small con- 
tribution due to equal-flavour quark scattering 
with t O u channel interference. The different 
gluon radiation pattern expected for Higgs bo- 
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Figure 1. Rapidity yj 3 of the third hardest jet, 
i.e. the one with highest p T after the two tagging 
jets. 



son production via VBF compared to its major 
backgrounds (ti production, QCD WW + 2 jet 
and QCD Z + 2 jet production) is at the core of 
the central-jet veto proposal, both for light 25 
and heavy [27] Higgs boson searches. A veto of 
any additional jet activity in the central-rapidity 
region is expected to suppress the backgrounds 
more than the signal, because the QCD back- 
grounds are characterized by quark or gluon ex- 
change in the {-channel. The exchanged partons, 
being colored, are expected to radiate off more 
gluons. 

For the analysis of the Higgs boson coupling to 
gauge bosons, Higgs boson + 2 jet production via 
gluon fusion may also be treated as a background 
to VBF. When the two jets are separated by a 
large rapidity interval, the scattering process is 
dominated by gluon exchange in the {-channel. 
Therefore, like for the QCD backgrounds, the 
bremsstrahlung radiation is expected to occur ev- 
erywhere in rapidity. An analogous difference in 
the gluon radiation pattern is expected in Z+2 jet 
production via VBF fusion versus QCD produc- 
tion [28 . In order to analyze this feature, we 
have plotted in fig. [1] the rapidity of the third 
jet, the one with highest p T after the two tagging 
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Figure 2. Comparison between the rapidity 
distribution of the third jet, yjz, obtained us- 
ing P0WHEG followed by the shower of HERWIG (in 
black), and HERWIG alone (in dark grey). The two 
curves have been normalised such that the total 
cross sections agree at NLO level. 



Figure 3. Comparison between the rapidity 
distribution of the third jet, y^, obtained us- 
ing P0WHEG followed by the shower of PYTHIA (in 
black), and PYTHIA alone (in dark grey). The two 
curves have been normalized such that the total 
cross sections agree at NLO level. 



jets. The conclusion that can be drawn from the 
figure is that the third jet in VBF prefers to be 
emitted close to one of the tagging jets, while, in 
gluon fusion, it is emitted anywhere in the rapid- 
ity region between the tagging jets. Thus, at least 
with regard to the hard radiation of a third jet, 
the analysis of refs. |29I30I31] is confirmed. The 
distributions, obtained using PDWHEG interfaced to 
HERWIG and PYTHIA, are very similar and turn out 
to be well modeled by the respective distributions 
of the NLO jet. 

For comparison, in figs. [2] and |3j we have plotted 
the same distribution as obtained with HERWIG 
and PYTHIA alone. Both the two plots are nor- 
malized such that the total cross section is correct 
at NLO. In fig. El it is clear thatthe third jet pro- 
duced by HERWIG underestimates the actual third- 
jet rapidity, even if the shape is similar. PYTHIA 
has instead a completely different behavior: here 
not only the normalization is wrong, but the pre- 
dicted 7/j3 distribution is not correct. This is of no 
surprise, since both HERWIG and PYTHIA produce 
(and resum) correctly only jets in the collinear 
region (collinear with respect to the two tagging 



jets), so that we do not expect any distribution 
involving the third jet to be correct. 

4. Z + 1 jet production 

In the early LHC phase, Z + 1 jet production 
is a promising process for jet calibration. Both 
Z and W production in association with jets are 
important sources of missing energy signals, and 
W in association with jets is an important back- 
ground to many new physics searches. The Z 
production process, with the Z decaying into lep- 
tons, gives a distinct signature that allows to 
perform stringent tests of the production mech- 
anism. These tests have been performed exten- 
sively at the Tevatron, both at CDF [32] and at 
DO [33l34|35l36j . They are carried out by cor- 
recting the measured quantities to the particle 
level, according to the recommendations devel- 
oped in the 2007 Les Houches workshop [37] , and 
then compared to NLO theoretical calculations, 
corrected for showering and underlying event ef- 
fects. These corrections, in turn, are extracted 
from SMC programs. 



lO" 



T3 



b 



10 1 



> 

u 
\ 

& 10° 



10" 



10" 



~~i 


LHC 14 TeV 


- \ 


Z 




Z+lj (p? orn > 5 GeV) 


- 


1 merged Z e Z+lj 

--I — i Mr - Mf — Pt 

i ^^ 


- 


'""--T^ -1 ^ 



100 



„ J1 
Pt 



200 
[GeV] 



300 



400 



> 

u 
s 

X! 



T3 



b 
T3 




Figure 4. The p T distribution of the first jet, 
in single Z production (blue dashed curve), in 
Z + 1 jet production (black solid curve) and of 
the merged sample (red dotted curve). 



Figure 5. The p T distribution of the second jet, 
in single Z production (blue dashed curve), in 
Z + 1 jet production (black solid curve) and of 
the merged sample (red dotted curve), when an 
angular separation Rjj > 3 is imposed between 
the two leading jets. 



Very recently we have implemented Z + 1 jet in 
the PDWHEG BOX. As a matter of fact, the P0WHEG 
BOX framework was developed with the Z + 1 jet 
process as its first example process. It is clear 
that, by using this tool, the comparison of the 
theoretical prediction with the experimental re- 
sults is eased considerably, and is also made more 
precise. Rather than estimating shower and un- 
derlying event corrections, using a parton shower 
program, these corrections can now be applied di- 
rectly to the hard process in question, yielding an 
output that can be compared to the experimen- 
tal results at the particle level. We have carry out 
this task and compared our results to the Teva- 
tron findings. We remark, however, that a fur- 
ther improvement to this study could be carried 
out, by using the P0WHEG program to generate the 
events that are fed into the detector simulation, 
and are directly compared to raw data. We are 
not, of course, in a position to perform such a 
task, that should instead be carried out by the 
experimental collaborations. 

Since a POWHEG implementation of Z produc- 
tion has already appeared in the literature [10j . 
a natural question arises, weather we can now, 
with the new Z + 1 jet calculation, build a merged 



sample of the two programs, capable of describ- 
ing Z production in the low, as well as in the 
high transverse-momentum region. We have pre- 
sented, for the first time, a minimalistic approach 
to such a merging, by putting together indepen- 
dent samples of Z and Z + 1 jet events in the 
most straightforward way, described in details in 
a forthcoming paper |16) . 

In fig. [4] we have plotted the p T distribution of 
the first jet, in single Z production (blue dashed 
curve), in Z -\- 1 jet production (black solid curve) 
and of the merged sample (red dotted curve). As 
can be seen in the figure, the merged distribu- 
tion smoothly interpolates between the single Z 
behavior at small p T , where the Sudakov form 
factor plays an important role in resumming the 
divergent collinear singularities, and the Z + 1 jet 
curve, at high p T , where the presence of the sec- 
ond hard jet in the Z + 1 jet sample, hardens the 
distribution of the first jet. Similar interpolating 
behavior is seen for the p T of the second jet, in 
fig- El Here, an angular separation Rjj > 3 is 
imposed between the two leading jets, in order 
to separate them further to better illustrate the 



behavior of the merged curve. 
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